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Abstract 

We construct a model in which the neutrino Dirac mass terms are of order 
the electron mass and the seesaw mechanism proceeds via right-handed neutrinos 
with masses of order TeV. In our model the spectra of the three light and of the 
three heavy neutrinos are closely related. Since the mixing between light and heavy 
neutrinos is small, the model predicts no effects in pp and pp colliders. Possible 
signatures of the model are the lepton-number-violating process e~e~ — t- H~H^, 
where H~ is a charged scalar particle, lepton-flavour-violating decays like — )• 
e~e~e~^, or a sizable contribution to the anomalous magnetic dipole moment of the 
muon. 
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1 Introduction 



The discovery of neutrino oscillations has established the existence of non-zero neutrino 
masses in the sub-eV rangejl] An appealing way of explaining such small neutrino masses 
is the seesaw mechanism in which gauge-singlet right-handed neutrino fields with 
Majorana mass terms are added to the Standard Model (SM). Those mass terms are 
not generated by the Higgs mechanism and may, therefore, be much larger than the 
electroweak scale. Let be the neutral members of the left-handed leptonic SM gauge 
doublets and Md and Mr the mass matrices of the fermionic bilinears urUl and unCuJ^, 
respectively. (C is the charge-conjugation matrix in Dirac space.) We denote the mass 
scales of Mjy and Mji by m^) and mji, respectively. If is much larger than m^, then 
an effective Majorana mass matrix 

= -MIM]^^Md (1) 

is generated for the z/^. According to ([1]), the scale my of A^j, is related to m/5 and m^j 
through niy ~ m|, / m^. The mixing between the light and heavy neutrinos is of order 
^d/itT'R, hence very small. 

We experimentally know that is in the range of 0.1 eV to leV, if 171,^ indicates 
the order of magnitude of the largest of the light-neutrino masses, but, in the framework 
of the seesaw mechanism, tud and rriR remain a mystery. Since is the neutrino 
counterpart of the charged-fermion mass matrices, mo may vary in between 100 GeV 
(which is both the electroweak scale and the top-quark mass scale) and 1 MeV (the scale 
of the electron mass and of the up- and down-quark masses), and this spans a range of 
five orders of magnitude. As tur ~ mjj / m^, to these five orders of magnitude in ttid 
correspond ten orders of magnitude in mR. If rriD ~ 100 GeV then ttir ~ 10^^ GeV; this 
is definitely below the typical GUT scale 10^^ GeV — identifying rriR with the GUT scale 
would make the neutrino masses too small. If instead ~ m,-, the mass of the tau 
lepton, then rriR ~ 10^ GeV. If one wants to incorporate leptogenesis [S] into the seesaw 
mechanism, then the appropriate rriR would rather be 10^^ GeV, which lies in between the 
two previous estimates. Finally, we may as well use mo ^ t^e, the electron mass, and then 
ruR ^ 1 TeV — this was noticed, for instance, in [6]H Such a low uir has the advantages 
that it coincides with the expected onset of physics beyond the SM and that it might 
produce testable effects of the seesaw mechanism at either present or future colliders. 

The possibility that m^) ~ rrie and ttir ~ 1 TeV is the starting point of this Letter. 
We construct a seesaw model in which the vacuum expectation value (VEV) responsible 
for the mass matrix Mo is of order rrie. In our model, which is inspired by ^ and [10], 
each charged lepton a {a = e, /i,r) has its own Higgs doublet (pa, whose VEV generates 
the mass rria- On the other hand, there is only one Higgs doublet 0o which has Yukawa 

"'^Cosmological arguments [1 and the negative result of the direct search for neutrino mass in tritium 
/? decay [2] are also crucial in eliminating the possibility that the neutrino masses may be higher than 
about one eV. For a review on neutrino masses see |3 . 

^Of course, the simplest way to obtain mu ^ rrie is to assume tiny Yukawa couplings, as was done for 
instance in [7] ; this is the opposite of what we do in this Letter — we discuss a scenario with Yukawa cou- 
plings of order 0.01-1. Another avenue which has been pursued are technicolor models with a suppressed 
rriD and toa, ^ 10'^ TeV [8^; however, in those models rajj is not claimed to be as low as mg. 
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couplings to the and is, therefore, responsible for M£,. We furthermore make use of a 
mechanism, first put forward in [11] and later extended in [12] , for suppressing the VEV of 
(pQ-. this doublet has a positive mass-squared /xq (in the scalar-potential term fio4'o4>o) and 
its VEV is triggered by a term 4>l4>o in the scalar potential. Since 0e is the Higgs doublet 
which gives mass to the electron, it must have a very small VEV, and this explains the 
smallness of the VEV of 0o- 

2 The model 

Multiplets: The gauge-S'?7(2) multiplets of our model are the following: 

• Left-handed lepton doublets D^a = {^La, cti)'^ and right-handed singlets i/rq, 
{a = e,/i,r). 

• Four Higgs doublets <^o and 

We shall use indices k, I running over the range 0, e, n, r. The VEV of the neutral 
component of (pk is denoted Vk- 

Symmetries: The family symmetries of our model are the following: 

• symmetry e ^ fi ^ r ^ e|f| 

• Three symmetries [9] 

Z„ : -DLa, OtR -> -Ur, Vrc, -l^Ra- (2) 

Notice that (pa does not change sign under z^. The symmetries may be inter- 
preted as discrete lepton numbers. 

• Three Z2 symmetries [9] 

Z^ : Or ^ -Or, (pa ^ (3) 

The z^ are meant to ensure that aR has no Yukawa couplings to the 0/? (/3 7^ a). 

Yukawa Lagrangian: The multiplets and symmetries of the theory lead to the Yukawa 
Lagrangian 

'Cy = -yi DLaaR(pa - ^2 XI ^Lal^Ra {ir2(p*o) + H.C. (4) 

a a 

This has, remarkably, only two Yukawa coupling constants. The Lagrangian (jl]) enjoys 
the accidental symmetry 

Z^ : 00 -00, T^Re -T^Re, T^Rp "I^R/., l^Rr -l^Rr- (5) 

■^If we go beyond our main aim of achieving a light seesaw scale and furthermore want to obtain the 
predictions of maximal atmospheric-neutrino mixing and a vanishing reactor-neutrino mixing angle, then 
one may extend this Z3 symmetry to the full S'3 permutation symmetry of e, and r |13'. 
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Charged-lepton masses: The mass of the charged lepton a is = \yiVa\- Therefore 
rrie '■ : 171^- = \ve\ '■ : \Vr\- In our model one cannot obtain a small rrig by just tuning 
some Yukawa couplings — one really needs a small VEV Vg. For instance, even if there 
were no further Higgs doublets in the theory and Vr by itself alone saturated the relation 

v'^J2\^kf<{l7AGeV)\ (6) 

k 

one would still need \yi\ ~ 0.01 because m,- ^ 1.78 GeV. We would then have \ve\ ~ 
50 MeV. On the other hand, we may assume that there are in the full theory some extra 
Higgs doublets other than the four (pk, which give mass to the quarks — in particular, 
for the top-quark mass a doublet with a large VEV is mandatory. Then \vr\ may be 
significantly smaller than 174 GeV and, accordingly, \ve\ will be significantly smaller than 
50 MeV too; in particular, |t>e| ~ me is possible. 



Soft symmetry breaking: We assume that the symmetry and the three Zq, sym- 
metries are softly broken in the dimension-three neutrino Majorana mass term^ 



(7) 



We furthermore assume that the symmetries are also softly broken, now by the 
dimension-two terms 0^0/ {k 7^ /) in the scalar potential. However, we shall assume 
that the combined symmetry 

Zg O Zj, : 00 -00, 4>e -0e, Gr -Cr, Ur^ -l^Re, l^Rf, 

is broken only spontaneously. Then the scalar potential is 

V = J2 /^fc0i0fc + ((j)i(j)k] 



•e ro4>e + fJ'fiT 0j,0r + H.C. 



(9) 



The quartic couplings in V, those with coefficients generically represented by the letter 
A, obey all the family symmetries of the model. Because of the couplings X'^i {(f)\(f)i^ , the 
only U{1) symmetry of V is the one associated with weak hypercharge; therefore, there 
are no Goldstone bosons in the model. 



Suppression of vq: If /iq is positive, then vq will induced by Ve and by the first term in 
the last line of ([9]): 

Vo ^ -Ve . (10) 

we want to predict maximal atmospheric mixing, then we must assume an S3 instead of a Z3 family 
symmetry and, furthermore, assume that the subgroup of S3, the /z-r interchange symmetry, is preserved 
in the soft breaking P1ITU1IT5]. 
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We envisage the possibility that I1/2I and, possibly, also \yi\ are of order 1, because that 
would enhance scalar effects and the experimental testability of our model, cf. sections 3 
and 4 below. Still, it is possible, as we mentioned earlier, that \yi\ ~ 0.01 and \ve\ ~ 
50MeV. However, even in that case |fo| could easily be much smaller than \ve\, as we 
pondered in [I2]. Indeed, |/ioe| could naturally [H] be small, and there is no reason why 
fio should not be rather large, maybe even of order TeV. Then |wo| might be much smaller 
than \ve\ — this is the mechanism that we envisaged in the introduction. Thus: 

• If there are in the theory some Higgs doublets beyond the four (pk, the Yukawa 
coupling constant yi in (jlj) may be of order one and then \ve\ ~ me- In that case, |/ioe| 
and /io may be allowed to be of the same order of magnitude and \vo\ ~ l^e] ~ me- 

• If there are only the four then yi is much smaller than one and \ve\ ~ 100 mg. 
In that case, we may naturally assume |/ioe| -C fio because the theory acquires the 
extra symmetry Zj, when /xoe = 0. We might then still obtain |fo| ~ ""^e- 



Lepton mixing: The neutrino Dirac mass matrix is in our model proportional to the 
unit matrix: 

Md = y;vol. (11) 
Therefore, the lepton mixing matrix U, which diagonalizes Aii,, also diagonalizes M^: 

U'^MuU = diag(mi,m2,m3) , (12) 
U^MrU* = -exp[2zarg(|/2*t;o)] diag(Mi,M2,M3), (13) 

the rrij and Mj {j = 1,2,3) being, respectively, the light- and heavy- neutrino masses. 
Consequently, there is a close relationship between the spectra of the light and heavy 
neutrinos: 

M, = (14) 
Following our rationale, we shall assume \y2Vo\ ~ 1^^. 



3 Possible collider effects 

In our model we assume ~ "^e ~ 1 MeV while mj? ~ 1 TeV. Otherwise our model is a 
normal seesaw model, therefore in it the mixing between the light and the heavy neutrinos 
is of order m£,/ mji ~ 10~^. This small mixing suppresses most possible signatures of heavy 
right-handed neutrinos that have been considered in the literature |T5]. For instance, the 
Drell-Yan production of a virtual W boson and its subsequent decay W^* — )• i^Nj, 
where £^ is a charged lepton and Nj a heavy Majorana neutrino, is negligible because it 
is suppressed by the mixing between the light and heavy neutrinos. 

Suppose that the Higgs doublet (pq (which may be one of our four doublets or an 
additional one) couples to the quarks and, in particular, generates the top-quark mass. 

^In that case we might envisage a scenario in which the VEVs of (f>r and 0^ would be responsible, 
respectively, for the masses of the up- type and down- type quarks. Other possibilities may of course also 
be considered. 
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Then we may envisage the Drell-Yan production of a virtual 
6q -)■ 00 and the decay 0o ~^ 



at the LHC, followed 



J' 

However, since the VEV of 0^ is necessarily large and the VEV of (pQ is very 



finally leading to heavy-neutrino 

;c 



will in general be small, unless we invoke finetuning in the scalar 



by the transition 
production 
small, the mixing 
potential. 

In contrast to what happens at hadron colliders, in an electron-electron collider the 
interesting lepton-number- violating process e~e~ — )• H~H~ might occur [6j. This process 
is due to the Majorana nature of the heavy neutrinos and is one of the processes, other 
than neutrinoless decay, via which it might be possible to probe lepton-number viola- 
tion p!6]. Let us suppose for simplification that H'^ = (pQ, the charged component of the 
scalar doublet ^q. Then, the relevant Yukawa Lagrangian for e~e~ — ?■ H~H~ is given by 



£ (e~e~ H-R-) = y^H^ ^ UI^N^Plc + H.c, 



(15) 



where Pl is the negative-chirality projection matrix. This leads to the total cross section 
(see [6] for a special case) 



a [e-e- H H 



(16) 



where s is the square of the energy of the e e system in its centre-of- momentum reference 



frame, /3 = (1 — 4m\j/ s)^^"^ 
function / is given by 



[niH is the mass of H ) and Wj = 1 — 2rn?fj/s + 2M|/s. The 



/ (a, h) 




a In 



fe + 1 



h ^ a, 
b = a. 



(17) 



Notice that the cross section f|T6|) depends on the Majorana phases of the products 
{UejU*f^)'^ {j 7^ k). In figure [1] we have plotted cr (e~e~ — t- H-R-) as a function of the 
light-neutrino mass mi in a number of cases. Following our rationale, we have taken 
\y2Vo\ = rrif, in equation f lT^ . On the other hand, in equation flTB]) we have taken \y2\ = 1, 
bearing in mind that a {e~e~ — )■ H~H~) depends very strongly on this Yukawa coupling. 
As for the U^j matrix elements, we have used the values |t/ei|^ = 0.7, \Ue2^ = 0.3 and 
Ue3 = 0, which almost coincide with the present best fit [T7]; this has the advantage 
that then the third-neutrino mass ma and the type of neutrino mass spectrum become 
irrelevant — we only have to take into account the experimental value of m| — m^, which we 
fixed at 8 X 10~^ eV^. We have moreover assumed that the phase of {UeiU*2)'^ is zero — this 
choice maximizes o" (e~e^ — ?■ H-H-). In figure [1] this cross section is given in units of 
o-QED = 47raV (3s), with a = 1/128. 

In the limit M| ^ s, mfj for all j = 1, 2, 3, one obtains 



2/3^ 

WjWk 



2MjMl ' 



:i8) 
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Therefore, in that hmit 



/3m 



ai^e-e- ^H-H-)^-—^, (19) 



where m^^ = J2j f^j (Uej) is the effective mass measured in neutrinoless decay. The 
approximation (fT9|) indicates a close relationship and, indeed, an approximate propor- 
tionality between cr(e~e~ H^H^) and m|^. Equation (fT9l) overestimates the cross 
section: a (e^e^ — t- H^H^) is much smaller than indicated by that approximation when- 
ever mi ^ 0.5 eV. On the other hand, for mi ^ 0.1 eV the approximation (IT^ becomes 
quite good, but at the same time the cross section becomes small. With the values of U 
used in figure [H for mi = 0.1 eV, y/s = 10^ GeV and mn = 120 GeV, equation f|T9|) gives 
a cross section about 14% too large; for smaller s or larger mn the discrepancy is smaller. 
Note that a cross section of about 10~^ x ctqed is still considered reasonable for detection 
of e~e~ — )■ H~H~ at an e~e~ collider [6]. 

4 Possible non-collider effects 

We next investigate whether large Yukawa couplings yi and ?/2 in equation (jl]) might 
induce measurable effects in non-collider physics. 
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4.1 The magnetic dipole moment of the muon 

A promising observable is a^, the anomalous magnetic moment of the muon. There is a 
puzzling 3a discrepancy [18] 

^exp _ ^SM _ 255(63) (49) x lO"" (20) 

between experiment and the SM prediction for that observable. In our model there are 
contributions to from one-loop diagrams involving either charged or neutral scalars. 
We consider the latter firstly. The fields 0^ written in terms of the physical neutral 
scalars (6 = 2, . . . , 8) as 
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VkbSh 



^l = v. + T.-^. (21) 



b-- 



t V2 ■ 



The scalar S\ corresponds to the Goldstone boson eaten by the gauge boson Z° and is, 
therefore, unphysical. The complex 4x8 matrix V is such that 



ReV 
ImV 



(22) 



is orthogonal^ therefore, iV^fel^ < 1. For b > 2, let rub be the mass of 5° and 5b = 
mj^/ml <^ 1. The contribution of the physical neutral scalars to a^^ is given by 

lAI /3 



E 

b=2 



b=2 

( X. [ \ A.\ /Q \ 1 

s:2 



167r2 



2 + In 5fe ) Re Ab 



+ 0i5t}, (24) 



where = l/iV^^. One sees that there is a term (5f,ln(5b /(IGtt^) which is ~ 10^^ when 
rrih ~ 100 GeV. Thus, if one assumes \ Af,\ ~ 1 then \a^ (0°)| might well be very large. 

We proceed to analyze the diagrams involving charged scalars in the loop. The three 
left-handed neutrinos and the three right-handed neutrinos are written in terms 
of the six physical Majorana neutrinos Xi (^ = 1; • • • ? 6) as 

6 6 

l^La = E Rc^iPhXh l^Ra = J2 ^aiPRXi, (25) 
i=l i=l 

where Pl,r are the chirality projectors in Dirac space and R, S are 3x6 matrices. The 
matrix 

(?) 



(26) 



is 6 X 6 unitary [20j. The fields (p^ are written in terms of the physical charged scalars 
S+ (a = 2,3,4) as 

0+ = 5:W,X- (27) 

a=l 



^For more details on this notation see 
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The scalar 5*^ corresponds to the Goldstone boson eaten by the gauge boson and is, 
therefore, unphysicaL The complex 4x4 matrix U is unitary [19j . 

Let rrii be the mass of the physical neutrino Xi ^^^y for > 2, let rria be the mass of 
S^. The contribution of the diagrams with charged scalars to is given by 



EE 



dx 



m. 



IQtt'^ ^2 i^i ■^^ m^x^ + {ml. — m?^ — mf^ x + m: 



X \^(\yiRfj,iV{^j,af + \y2Sf,iUoaf) [x^ - 



X 



X — X 



:28) 



The three light neutrinos Xi,2,3 have masses much smaller than m^, hence these masses 
are negligible. For those neutrinos the matrix elements S^j^i ~ 10^^ are also negligible. 
One then has 



light neutrinos 



967r' 



E H 



2 



(29) 



a=2 



where 6a = rn? Jrn?^. For rUa ~ 100 GeV one thus has 



-10-9 



2 

11 ) 



"/!/ ■■"a' ^ ■■"a . „„„ yyfy yiigjit neutrinos 

which is not very significant. The three heavy neutrinos X4,5,6 have masses comparable to 
those of the charged scalars. For those neutrinos the R^i ~ 10"^. One then has 



heavy neutrinos 



/ \y2S ^ihlQa\^ 1 

Re [ylylR,,S^dA;M^ 




where 



X 



+ 



X 



a; In X 



3(x-l) 2(x-l)' {x-lf {x-lf 



1 

-3 + 



" 1 1 

mx\ X 

6 



(30) 

(31) 
(32) 



For Yukawa couplings of order 1, and for heavy neutrinos with masses of order 1 TeV, 
one obtains (0^)hcavy neutrinos ~ 10^^°, which is smaller than the experimental error in 
equation (|20l) and hence irrelevant. 

So one concludes that the largest non-standard contribution to in our model is in 
principle (0°), which is proportional to \yi\^ and may be as large as 10"^ if \yi\^ ~ 1. 
One must invoke either a small Yukawa coupling yi or cancellations between the contri- 
butions of the various neutral scalars for our model not to give too large a contribution 
to a^. One may also view fl21|) as a possible way to explain the discrepancy fl2U]) . 
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4.2 The decay fi 



e e ' e 



In our model, there are also lepton-flavour-changing processes due to the soft breaking of 
the family lepton numbers [20] . The process most likely to be observable is in general [2^ 
//~ — )■ e~e^e~^ which is mediated by /i~ — )■ e~S'°* if the Yukawa couplings are large. 
Using the formulas in ^20j we obtain the approximate upper bound 



br(, 



/i e e e 



< 10 



-5 I |4 

\yiy2\ 



. b 



2 (100 GeV)^ 



mt 



where 



3 

i=l 



(33) 



(34) 



fl being an arbitrary mass parameter on which F finally does not depend. Clearly, if 
yi, 1/2 and F are all of order one, then BR — > e~e+e~) will be much too large in our 
model, since experimentally that branching ratio cannot be larger than 1.0 x 10~^^ at 90% 
confidence level [18]. One possibility to avoid this problem is by making \yi\ ~ 0.01, which 
is possible as we have seen in section [21 and would have the further advantage of also sup- 
pressing (0°), while keeping ?/2 at order one in order not to suppress a (e~e~ — t- H^H^). 
Another possibility is to assume that the mass spectrum of the light neutrinos is almost 
degenerate, which in turn renders the mass spectrum of the heavy neutrinos almost de- 
generate too. With this condition and assuming [/gs to be exactly zero, we estimate 



mf 



(35) 



Taking, for instance, mi = 0.3 eV, one obtains \F\ ~ 10"'^, which is sufficient to make 
equation (l33ll compatible with the experimental bound. 

Therefore, our model may be compatible with the present experimental bound on 
BR(/i~ — )■ e~e~^e~), if the mass spectrum of the light neutrinos is almost degenerate or 
the Yukawa coupling yi is of order 0.01. On the other hand, the Yukawa coupling y2 may 
perfectly well be of order one. 



5 Conclusions 

We have started in this Letter with the trivial observations that, in the seesaw mechanism, 
a mass scale of 1 MeV in the neutrino Dirac mass matrix corresponds to right-handed 
neutrinos in the TeV range, and that the scale 1 MeV might be provided by the electron 
mass. We have then constructed a simple model, which extends the SM with three right- 
handed neutrino singlets and four Higgs doublets, in which that observation is realized. 
Our model has the following properties: 

1. Each charged lepton a has a corresponding Higgs doublet (pa which, through its 
VEV Va, generates the charged-lepton mass = 

^The process /x~ — >■ e~7 usually is more suppressed than jj," — >■ e~e+e^ in this type of models |20) . 
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2. The Dirac mass matrix M/) of the neutrinos is generated by another Higgs doublet, 
00, whose VEV Vq is induced by the VEV such that Vq ~ Ve or smaller, hence 

Vq oc rrif,. 

3. In the appropriate weak basis, the mass matrix of the charged leptons is diagonal 
while Mo is proportional to the unit matrix; this yields the simple relation Mj oc 

/ nij between the masses Mj of the heavy Majorana neutrinos and the masses 
nij of the light Majorana neutrinos. 

4. Moreover, the diagonalization matrix of the light-neutrino mass matrix — which is 
just the lepton mixing matrix — and the diagonalization matrix of the heavy-neutrino 
mass matrix are complex conjugate of each other. 

5. The mixing between light and heavy neutrinos is small, of order 10~^. 

The last property prevents heavy-neutrino production in pp and pp colliders. However, at 
an e~e~ collider one might test the mechanism of the model by searching for the process 
e~e~ — i- H~H~, where H~ is a charged scalar; this process is somehow a high-energy 
analogue of neutrinoless decay. The charged scalar would mainly decay into a light 
neutrino and the electron, if the heavy neutrinos are heavier than the scalar. Thus the 
signal would amount to e~e~ plus missing momentum. At an e~^e~ collider one should, 
of course, look for e+e^ H^H~. 

If the Yukawa coupling yi of our model is of order one, then the anomalous magnetic 
moment of the muon may be too large, and the branching ratio of the flavour- changing 
process ^~ — )■ e~e~^e~ as well; in this case the suppression mechanisms discussed in 
section H] have to be invoked. The simplest way to make our model compatible with 
experimental constraints is by choosing \yi\ ~ 0.01. However, this does not impede the 
process e~e~ — t- H~H~ , which proceeds through a different Yukawa coupling ?/2- 
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